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Abstract 

A non-gravitational, non-linear interaction between dark matter and dark energy may result 
in a future evolution of the Universe which differs from that of the standard ACDM model. In 
particular, the ratio of the energy densities of dark matter and dark energy may approach a stable 
finite value. For a special case we find a corresponding analytic solution for the interacting two- 
component dynamics which is consistent with the supernova type la (SNIa) data from the Union2 
set. For a broader class of interactions without analytic solutions, a dynamical system analysis 
classifies stationary points with emphasis on their potential relevance for the coincidence problem. 
Asymptotically stationary solutions of this kind require a phantom-type "bare" equation of state 
of the dark energy which, however, does not lead to a big-rip singularity. 
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I. INTRODUCTION 



Investigating the properties of the cosmological dark sector has become one of the major 
activities in physics since the detection of the accelerated expansion of the Universe more 
than a decade ago [1]. According to the most accepted interpretation, based on Einstein's 
General Relativity, our Universe is dynamically dominated by two so far unknown compo- 
nents, dark matter (DM) and dark energy (DE). The latter contributes roughly 72% to the 
total energy budget, the former about 23%. Only about 5% are in the form of conventional, 
baryonic matter. DE, a substance equipped with a sufficiently large negative pressure ac- 
counts for the accelerated expansion, DM is needed for successful cosmic structure formation. 
Alternative and complementary support for this interpretation comes from the anisotropy 
spectrum of the cosmic microwave background radiation [2J , from large-scale-structure data 
[3], from the integrated Sachs- Wolfe effect [1], from baryonic acoustic oscillations [5] and 
from gravitational lensing [6j. The preferred model is the ACDM model which also plays 
the role of a reference model for alternative approaches to the DE problem. Because of the 
cosmological constant problem in its different facets and the coincidence problem, i.e. the 
question, why the ratio of the energy densities of DM and DE is of the order of unity at the 
present epoch, a host of alternative models has been developed in which the cosmological 
term is dynamized. Overviews of the situation can be found, e.g., in |71-E]- While in the 
ACDM model and in most of the alternative approaches DM and DE are considered as 
independent components of the cosmic medium, there exists a line of research that admits 
a coupling between both dark components. Such coupling does not contradict the overall 
energy-momentum conservation. Interacting models of this type give rise to a richer dynam- 
ics and are particularly useful to address the coincidence problem. Models, applicable to an 
interaction between DE and DM were introduced by Wetterich [10]. Meanwhile there exists 
a still growing body of literature on the subject -see, e.g., [TTH2n] and references therein. A 
general problem here is the choice of the interaction term. Since neither the physical nature 
of DE nor the physical nature of DM are known, there does not exist a sound microphys- 
ical motivation for a specific interaction either. Therefore, approaches to interacting dark 
energy are largely phenomenological. Even though, it has been argued that, being unaware 
of the possibility of an interaction between both dark components, may result in a misled 
interpretation of observational data [12]. Most of the interactions studied so far are linear 
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in the sense that the interaction term in the individual energy balances of the components 
is proportional either to the DM density or to the DE density or to a linear combination of 
both densities (for a recent analysis see, e.g., pR]). Since systems with interactions admit 
analytical solutions only in special cases, several authors resorted to a dynamical system 
analysis (see, e.g., [301132] for the general background) to obtain a qualitative picture of the 
long-time behavior of the cosmological dynamics |TiH20] . Occasionally, non-linear couplings 
were studied as well |2Tl |22]. A product coupling, i.e., an interaction proportional to the 
product of DM density and DE density was shown to be favored observationally over linear 
choices in [23]. Also from a physical point of view such type of coupling seems preferred. 
An interaction between two components should depend on the product of the abundances 
of the individual components, as, e.g., in chemical reactions. A dynamical system analysis 
for a specific non-linear interaction was performed and contrasted with observational results 
in [211 [25]. 

It is the purpose of this paper to investigate the long-time behavior for a simple two- 
component model with a number of non-linear interactions on the basis of Einstein's theory. 
We show that for a particular case there exists an analytic solution with a non-vanishing, 
positive limit for the ratio of the energy densities of DM and DE. In this case, the equation- 
of-state (EoS) parameter w of the DE is necessarily of the phantom type, i.e., w < — 1. 
Similar stationary points for a larger class of interactions are found with the help of a 
dynamical system analysis. We classify these points according to their possible relevance for 
the coincidence problem. The existence of these points requires w < —1 as well. However, 
while non- interacting models with constant w < —1 necessarily approach a singularity after a 
finite time [33l |3l] , the (non- linear) interaction quite generally prevents the cosmic evolution 
from a final big-rip. Moreover, independent of the details of the interaction, a positive total 
energy density in the critical points necessarily implies an energy transfer from DE to DM. 
We discuss attractors, stable focuses and centers as potential final states of the cosmic 
dynamics. The center-solution discussed in [25] is recovered as a particular case of our work. 

The paper is organized as follows. In section[TT]we present the basic relations of the general 
interacting two-component model. We clarify that the relevant critical points require a 
phantom-type EoS parameter and we show that a positive critical density is only compatible 
with a production of DM at the expense of DE but not with a transfer in the opposite 



direction. A broad class of non-linear interactions is introduced in section III , where we also 
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test a particular analytic solution against SNIa data from the Union2 set |35] . In section IV 



we perform a dynamical system analysis. A discussion of how the long-time limit might fit 



into a viable cosmological scenario is given in section |Vj Section |VT] summarizes the results 
of the paper. 

II. GENERAL ANALYSIS 

The dynamics of our present Universe is assumed to be dominated by DE and DM. The 
relevant field equations for the spatially flat, homogeneous and isotropic case then are the 
Priedmann equation 

3i/2 = 87rG (p^ + p,) , (1) 

and 

H = -AtxG {Pm + Px+Px) ■ (2) 

Here, pm is the energy density of pressureless DM and px is the density of a DE component 
with a pressure px- Throughout this paper we use units with c = 1. We assume that both 
components do not conserve separately but interact with each other such that the balance 
equations take the forms 

pm + 3ifpm = Q (3) 

and 

Px + ?>H{1 + w)px = -Q , (4) 

where w = is the EoS parameter of the dark energy. The sum of (j3j) and Q results in 
the total energy conservation equation 

p + 3H{p + p)=0, (5) 



where the total pressure equals the dark energy pressure, p = px- To address the coincidence 
problem it is convenient to in1 
characterized by the dynamics 



problem it is convenient to introduce the ratio r = ^ of the energy densities, which is 



(6) 



Pm Px 
Pm px _ 

Furthermore, it is also convenient to introduce an effective pressure quantity 11 by Q = 
— 3i7n and to replace the derivatives with respect to the cosmic time by derivatives with 
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respect to In a^, denoted by a prime, i.e. p = p'3H. Then the dynamics of the two-component 
system is given by 



Pm 
Pra 



-1 



n 

Pm 



Px 



[l + w) + 



n 

Px 



or, alternatively, by 



and 



w 



:i+r) 



w 



n 



rp 



(7) 
(8) 
(9) 



All the details of the interaction are encoded in the function 11. In the interaction-free limit 
n = 0, the stationary point = together with w = —1 corresponds to the de-Sitter space 
as the long-time limit of the ACDM model. 

The relevant critical points of Eq. ([8| are given by 
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w 



(10) 



where the subscript c denotes the critical point. Consequently, for positive values of r, the 
existence of a critical point requires an EoS parameter w < —1, i.e., DE of phantom type. 
This conclusion does not depend on the interaction. A non-zero stationary value for the 
ratio r can be interpreted as an alleviation of the coincidence problem. The condition r' = 
together with ^ and (10) provides us with 

w 



Pc 



l + w 



(11) 



In general, = Yic{Pci^c)- Therefore (11) is not an explicit relation for p^- Moreover, 
pc remains undetermined for a linear dependence of 11 on p. As will be shown in section 



IV below in more detail, this case is degenerate and does not admit a dynamical system 



analysis. On the other hand, for 11 oc p equation (|9]) decouples which will allow us to obtain 
analytic solutions of the system ^ and 

Since w < —1, a positive stationary energy density pc in (11) requires Tic < 0, equivalent 
to Qc > 0. Independent of the specific interaction (excluding only a linear dependence 
n oc p), the existence of the critical points and pc requires a transfer from DE to DM. 
We disregard here the critical points with Vc = Pc = and the unphysical Tc = —1 (pc = 0). 
We emphasize that as long as 11 oc p is excluded, the results for the critical points so far do 
not depend on the structure of 11. In the following section we shall consider a specific class 
of interactions. 
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III. A CLASS OF NON-LINEAR INTERACTIONS 



A. Structure of the coupling term 

As already mentioned, lacking a microphysically motivated interaction, one has to resort 
to phenomenological models. Notice that the knowledge of a reliable microphysical inter- 
action would be equivalent to already knowing the physical nature of DE. Our principal 
interest in this paper are non-linear interactions, i.e. interactions for which the effective 
pressure 11 in general is a non-linear function of the energy-densities of the components 
and/or the total energy density. Motivated by the structure 

r 1 
Pm = z—, — p and = — — p (12) 
1 + r 1 + r 

of the components, we consider the ansatz 

n = -7p-r" (1 + rY = -ip-^+'p-^pl- , (13) 

where 7 is a positive coupling constant with a dimension of p^^™. The powers m, n and s 
specify the interaction. For fixed values m, n and s the only free parameter is 7. A linear 
dependence of 11 on p corresponds to m = 1. While this case is not accessible to a dynamical 



system analysis (see the comments following eq. (11) and section IV below) it is perfectly 
admissible in the general dynamics (|8| and (|9]). Moreover, as already mentioned, it is exactly 
this case which will provide us with analytic solutions of the system ^ and ([9]). 

The effective interaction pressure 11 is proportional to powers of products of the densities 
of the components for the special cases s = —m, but we shall admit s and m to be arbitrary 
for the moment. Notice that every power of the total energy density p in the interaction 
pressure itself corresponds, via Friedmann's equation, to the square of the Hubble parameter 
(p oc H"^). This implies that the interaction quantity Q is not necessarily linear in the Hubble 



rate. For s = —m the ansatz (13) is equivalent to 



Q = 3H^p^-yi = 3H^p^r^ . (14) 



The ansatz ( 13 ) contains a large variety of interactions that have been studied in the liter- 
ature as special cases. This comprises, e.g., the models in [T514T71 [26H29] . 

As mentioned before, in most of the interacting models in the literature, the interactions H 
are assumed to be linear in either p^ or p^. The corresponding source terms are Q = 3'jHpm 



or Q = respectively, or a combination of both (see, e.g., [T5HT71 |29])- In our 

setting, the case Q = 3'yHpm is recovered for {m,n, s) = (1, 1, —1), while the combination 
(m, n, s) = (1, 0, —1) reproduces Q = 3'yHpx. 

In the following subsection we consider particular combinations of the parameters (m, n, s) 
which give rise to analytically solvable models with non-linear interaction terms. 



B. Analytically solvable models 



1. The case Q = 3H^ 



pmpx 



This example for an analytically solvable non-linear interaction model, covered by the 
ansatz (13), follows for (m, n, s) = (1, 1, —2). In such a case, equation ^ reduces to 



r' = r[w + ■y] 



which results in a power-law solution 



roa 



3(ui+7) 



(15) 



(16) 



The ratio r decreases for w + j < 0. Introducing (|16|) into (|8|), we find the energy density 



P = Poa 



1 + TQa 



3(«)+7) 



(17) 



1 + ro 

This case coincides with the interacting model studied in relying on an ansatz 

r = roa~^ for the energy-density ratio. This ansatz was proposed in [33] in order to address 
the coincidence problem. The correspondence is 7 = — (w + |) . For a <^ 1 we have p oc a~^, 
for a ^ 1 the behavior is p oc a~^^^~^^\ The ACDM model is recovered for ^ = 3 and w = —1, 
equivalent to 7 = 0. The densities of the components are 

' 1 + roa^^'"+^y 



Pr. 



PmOO- 



-3(1-7) 



and 



where p. 



Px — Pxao- 



-3(l+-t«) 



1 + ro 

1 + roa3("'+T) 
1 + ro 



7 
w-\-'Y 



(18) 



(19) 



mO 



'"0 



Po and pxo = j^Po, respectively. The non-interacting limit is correctly 

isity rati 

3(l+Tff) 



1+ro' " ' J-+''0 

reproduced for 7 = 0. For w + 7 = the energy- density ratio r is constant and 



r = To 



p = Poa 



(20) 



The model based on (17) has been analyzed in some detail in the literature [38| H0| WT\ . 
It represents a testable alternative to the ACDM model. Although the latter is largely 
consistent with observations, the data leave sufficient room for deviations from either ^ = 3 
or w = —1 which would correspond to a non- vanishing interaction, i.e., 7 7^ 0. 



2. The case Q = J.H-f' 



This case corresponds to a choice (m, n, s) = (1, 2, —2). It has the analytic solutions 



w 



and 



r = ro 



_3( I 

p = poa ^ 



{w + 7ro) a — "jro 

{w + 7ro) a~^^ + tq (w — 7) 



w{l + ro) 



The high-redshift limit of (21) is 



r — 7- 



\w\ 
7 



(a < 1) 



(21) 



(22) 



(23) 



For a <C 1, i.e. in the past, the ratio r becomes constant. In the opposite limit a 3> 1 on the 
other hand, we find r oc as in the ACDM case. Assuming w = —1, the energy density 
for small values of the scale factor behaves as 



p (X a i+T 



(a < 1) . 



(24) 



The interaction constant modifies the typical behavior at high redshifts. In the opposite 
limit a ^ 1, the energy density behaves as p oc a^^^^^"^^ which coincides with the corre- 
sponding dependence of the wCDM model. In this case, the interaction does not lead to a 
different future evolution of the Universe. 



3. The case Q = 3H^ 



The third and most interesting analytical solution corresponds to the choice (m, n, s) 
(1,0, —2). For w <0, i.e. w = —\w\, the solutions are 



\w\ I \w\ 



(25) 
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and 



P = PoO- 



-3 1- 



+ 7 + {\w\rQ - 7) a' 



-3\w\ 



\w\ (1 + ro) 

For 7 > 0, positivity of both r and p is guaranteed for 7 < \w\rQ: 



l™l+7 



7 < \w\rQ 



r > and p > . 



(26) 



(27) 



The ratio (25) scales as r oc a ^'"'1 for a <^ 1. For w = —1 this coincides with the scahng of 



its ACDM counterpart. In the far-future hmit, however, we have 

(a > 1) , 



r — )■ 



(28) 



i.e., the energy-density ratio remains finite, whereas it tends to zero in the ACDM modeL 
The energy density scales as a^'^ for a ^ 1, i.e., we recover an early matter dominated 
period. In the limit a ^ 1 one has 



p (X a 



-3 1- 



|iti|+7 



(a > 1) 



which generally does not correspond to a de Sitter phase. 



The solution (26) has the interesting special case 1 — 



|ui|+7 



(29) 



in which p tends to a 



constant for a ^ 1. Under this condition we have 



\w\ = \w\ + 



7 = \w\ {\w\ — 1) 



\w\ 



\w\ — 1 



The interaction constant 7 is directly related to the deviation of w from w 



energy density (26) may be written as 

\w\ + [ro - {\w\ - l)]a-3l'"l 



P = Po 



1 + ro 



1 



For r we find 



The limiting values for a ^ 1 are 



Poo = Po 



\w\ 



\w I 



and 



1 + ro 



Ityl — 1 



(30) 

■1. Then, the 

(31) 
(32) 

(33) 



(34) 



The dynamics results in stationary values for p and r. Notice that p > Poo and r > Too- 



Moreover, the limiting value ( 34 ) for r coincides with the stationary value ( 10 ) . In particular, 
we have again w < —1. To the best of our knowledge, this solution has not been considered 
before. The role of the interaction in the limiting cases is as follows. In the distant past, 
e.g. at high redshift, we have r ^ 1 and p ~ pm- Then, ~" ~ ~" oc ^ 1, i.e., the 
interaction is negligible and a matter dominated phase is correctly reproduced. On the 

7 and 



|n| _ |n| 



other hand, in the far future, r ^ 1 and p ~ p^ are valid. Consequently, ^ ~ ^ ~ 



M ~ 7 
Pm r 



. In this limit the interaction is crucial. In other words, the interaction is switched 



on during the cosmic evolution. In the following we check whether the solutions (31) and 



(32) with the final stationary values (33) and (34), respectively, are consistent with current 



SNIa observations. The crucial quantity is the Hubble rate that corresponds to the energy 



density (26): 



H = Hoa 



|uj|+7 



+ 7 + {\w\ro - 7) a 
\w\ (1 + ro) 

The deceleration parameter changes from q = \ for a ^ 1 to 



2 |tii|+7 



1 - 



3 w 


2 


\w 


+ 7 



(a > 1) 



(35) 



(36) 



Notice that for the case (30) the limiting value is g = —1, although Too > according to 



(34), i.e., different from the ACDM model there remains a non- vanishing matter fraction. 



For our statistical analysis we do not specify beforehand to the solutions (31) and (32) 



Our aim is to clarify whether the parameter combination (|30|) that corresponds to these 

^0 



solutions has observational support. The free parameters are \w\^ 7 and Vt^ 



1+ro ■ 



We performed a Bayesian statistical analysis on the basis of the SNIa data from the 
Union2 set [35], using the marginalization method for Hq developed in [12]. The results 
of this analysis are shown in Fig. [!} The best-fit parameter values are (7, w, ^o, x^mm) = 
(0.36l?;^^,-1.20l[^:f,0.39l[5:?^, 540.863). These results demonstrate, that within the la re- 



gion the solutions (31) and (32) with the final stationary values (33) and (34), respectively, 
are indeed consistent with the observational data. This may indicate a phenomenological 
solution of the coincidence problem with the help of non-linear interactions in the dark 
sector. 

For an arbitrary combination of the parameters m, n and s, analytical solutions of the 
non-linear system are hardly available. To get insight into the behavior of the system under 
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(a) 




-0.5 0.0 0.5 1.0 1.5 



y 



(b) 



Q - 3H 7 px"'p 




-1.8 -1.6 -1.4 -1.2 -1.0 -0.8 -0.6 
w 

(c) 

FIG. 1: Interaction model with (m, n, s) = (1, 0, —2). The contour plots denote the la, 2a and 3a 
regions, based on the Union2 data set. 



more general conditions we shall resort to a dynamical system analysis in the following 
section. 



IV. A DYNAMICAL SYSTEM ANALYSIS 

In this section we consider the dynamics in the vicinity of the critical points with the help 
of a dynamical system analysis. This analysis is based on the circumstance that, close to 
the critical points, the (generally unknown) solution of the non-linear system behaves as the 
solution of the system, linearized around the critical points (Hartmann's theorem and, for 
purely imaginary eigenvalues, the Center Manifold Theorem (see, e.g., [32j and (SO])- Using 
standard techniques (see, e.g., [311 [32]), the general characteristic equation for the critical 



11 



points is 



2 + w — w (1 + w) 



IT 



A + (1 + w + wdpU) = 



(37) 



Here 9^11 and dpU denote the partial derivatives of 11 with respect to r and to p, respectively. 



Eq. (37) has the solutions 



w{l + w) — - {2 + w) 



± 



2 + w — w {1 + w) 



IT 



-4(1 + 1;; + wdpU) 



(38) 



where we have to require 1 + w + wdpU ^ 0. In case these solutions are non-degenerate 
and real, they describe an attractor for \± < 0, an unstable critical point for A-t > and 
a saddle if A+ and A_ have different signs [32] • For complex eigenvalues X± = a + if3, it is 
the sign of a that determines the character of the stationary point. For a = the critical 
point is a center, for a < it is a stable focus and for a > it is an unstable focus. With 



the ansatz (13) for the interaction, the eigenvalues (38) are 



\± = l[2 + s + {1 + n + s) w] T ^(2 + s + (1 + n + s) w)^ + 4 (m - 1) (1 + w) 



(39) 

Since the analysis is valid only for non-zero eigenvalues, all cases with m = 1, corresponding 
to a linear dependence of 11 on p, are not covered by the classification mentioned before. 
Notice that both the linear interactions and the analytically solvable non-linear cases of the 
previous section have m = 1. For m = 1 equation ^ is decoupled from p and can be solved 
separately. 

For m 7^ 1 the general classification provides us with the following set of critical points: 



Attractor for m > 1 and s < 



Unstable for m > 1 and s > 



2+(l+n)w 



1+w 



-2 



2+(l+n)w) 
1+w 



+ 2 



l—m 
1+w 



l—m 
1+w 



Saddle for m < 1, for all n and s 



Center for m > 1 and 2 + s + (1 + n + s) w = for n > 1 



Stable focus for m > 1 and 2 + s + (l + n + s)w>0 



Unstable focus for m > 1 and 2 + s + (l + n + s)w<0 
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For the critical density we find (m 7^ 1) 
The corresponding values for the components are 

Pmc=[l\wY^'^{\w\-ir-'^]-- 



(40) 



and 



consistent with 



(41) 



(42) 



— 1 in (10). Consequently, the fixed points {pc,rc) are given 



by (40) and (10). For m > 1 the critical densities are proportional to a negative power of 



the interaction constant, i.e., the diverge in the limit 7 — i- 0. This limit corresponds to the 
big-rip singularity of dark energy models with constant EoS parameters w < — 1 [331 [3l]. In 
other words, any non-vanishing interaction of the type considered here is big-rip avoiding. 

In tables |T| |Tl] and III we present potentially interesting examples for different critical 
points, together with the allowed ranges for the EoS parameters. Fig. ^ shows the corre- 
sponding phase-space portraits. 

For a physical interpretation it is instructive to check the analytic solutions of the lin- 



earized system that underly (38) and (39) explicitly for special cases. 



(i) An attractor with {m,n,s) = (2,0,-2). This corresponds to 11 = — 7p^ (1 -|- r) 



The general expression (40) for the critical density reduces to 

\w\ {\w\ — 1) 



Pc 



1 



(43) 



Eq. ([9]) simplifies to 



r' = rw + •yp . (44) 
Let us introduce quantities / and g which describe small deviations from the critical point: 



P = Pc + f , r = r^ + g . 



(45) 



While in zeroth order r^w -|- = is valid, we have from (44) and ([8j), up to linear order 
in / and g, 



g' = gw + -ff , and /' 



1 
7 



{\w\-l) . 



(46) 
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m 


n 


S 


n 




W 


3 

9 





-1 


-7/j3/2(i + r)-i 




-1.5 <w<-l 


3 
2 


1 

2 


3 
2 


-7p3/Vl/2 (1 + ^)-3/2 




-1.125 < w < -1 


3 
2 


1 
2 


-1 


_^p3/2^1/2 ^yl 


^H^y/pPxPm 


-1.101 <w <-l 


2 


1 
2 


-3 




SHjPx^pPm 


-1.0625 <w <-l 



TABLE I: Examples for an attractor as critical point. 



Eliminating / yields 



g"+\w\g' + i\w\-l)g = 0. 



(47) 



This second-order equation with constant coefficients has the solution g = Qq exp [Ax] , where 
go is some initial value and A is determined by A^ + \w\\ + {\w\ — 1) = 0, i.e., Ai,2 = 
— ^ ± ^J~^^^^^-^-~i]w\~-^l) . Both solution are negative and, consistently, represent a special 
case of (39). Consequently, g decays exponentially with x. Recalling that x = Ina^, we have 



9 = 9oa 



3A 



and 



3A 

— ^oa 

7 



(48) 



Since A is negative, g and / decay with a power of the scale factor. It is the scale-factor 
dependence of (48) that is behind the curves in Fig. |2](a). 

The expression in ([8| represents the total effective equation of state of the cosmic 



medium. Close to the critical value 



l + Tc 



W 



—1 it is given by 

,3A 



■1 + 



9 



\w\ 



1 + r 

The deceleration parameter is related to by 

H 



1 + 



\w\ 



(49) 



-1 - 



if2 



g=-(l + 3^ 
^ 2 V 1 + r 



(50) 



In the critical point itself gc = ~1 is valid. In the vicinity of the critical point we have 



with g from (|48|). Both the total equation of state and the deceleration parameter q 



approach —1 with the power 3A of the scale factor. 
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FIG. 2: Phase portraits: in (a) the interaction Q = 37il;^^^= corresponds to an attractor for 
w = —1.08 and 7 = 1. The same interaction results in a stable focus for w = —1.2 and 7 = 1 in 
(b). In (c) the interaction Q = SjHp^px describes a center with w = —1.17 and 7 = 1. 
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TABLE II: Examples for a stable focus as critical point. 

(ii) A similar analysis can be made for for the dynamics around a center with {n, s) — 
(1,2), corresponding to 

n = -7p"^-VmPx ^ Q^SH^p'^-^PmPx . (52) 
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TABLE III: Examples for a center as critical point. Notice that the second case here corresponds 
to the model studied in [25l. 



In this case, the critical density is 



Pc 



\w\ 

7 



1 

m— 1 



(53) 



With the ansatz (45) we have, up to linear order, 

p--i = p--i l + (m-l 
The linearized system becomes 

—7 (m — 1) (1 + w) 



pc 



(54) 



9 



\w\ 
7 



m-2 
m — 1 



/ 



and 



w 



\w\ 
7 



(55) 



This is equivalent to the second-order equation /" — (m — 1) (1 + w) / = which, for 
m > 1, describes oscillations with a frequency oo = ^(1 — m) (1 + w), i.e., / = 
/o cos ^a/ (1 — m) (1 + w) where x = Ina^. For g we have an identical equation with 
solution g = (^osin — m) (1 + w) , i.e., the dynamics is described by ellipses in the 

f — g plane. This behavior is visualized in Fig. ^c) The change of the total energy density 
in ^ is given by 

- = - = -j^- (56) 

p p l"*^! 

For g > 0, i.e., if the ratio r is larger than the critical value, equivalent to an enlarged matter 
contribution, we have p' < 0, i.e., the total energy density decreases since the effective EoS 
parameter is larger than —1. On the other hand, for g < 0, there is an excess of dark energy 
and the dynamics is given by p' > 0. The rate of change of p is oscillating about the critical 
point. The period ■^/(l — m) (1 + w)x = 2ti corresponds to ap = exp 
m = 2 and w = —1.1, e.g., we find ~ 750 numerically. 



2tt 



. For 
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In the vicinity of the critical point we have 

With an oscillatory solution for g the deceleration parameter oscillates around g = — 1. This 
corresponds to an oscillation of the effective total equation-of-state parameter around 
— 1 as well. The cosmic medium as a whole oscillates between phantom- and non-phantom 
behavior, the phantom divide is crossed periodically. 

(iii) A mixed case is the stable focus, e.g., the case (m,n, s) = (|, ^, — l). By a similar 
analysis we find that the perturbations about the critical point behave as damped oscillations 



with a frequency (with respect to x) of 2 {\w\ — 1) — (l + and a damping factor 
^-3(i+u)/2)/2_ rjj_^g critical values —1 for the effective equation of state and the deceleration 
parameter are approached by a corresponding spiralling-in dynamics which is visualized in 
Fig.gb). 

(iv) The case m = n = s = Oisan example for a saddle. Here we have 11 = const. 



Consistent with the corresponding special case of (39), only one of the solutions approaches 
q = —1 with a power of the scale factor, the second solution is unstable. 

It is interesting to realize that the same interaction may result in different critical points 
for different ranges of the EoS parameter. The third example in Table [T] corresponds the 
same expression for Q as the fourth example in Table [TTj For the range —1.001 < w < —1 
the critical point is an attractor while it is a stable focus for —2 < w < —1.101. A similar 
situation occurs for the second case in Table [T] and the third case in Table [Tll 

V. AN EARLY MATTER-DOMINATED PHASE 

The focus in this paper is on the late-time behavior of the cosmological dynamics. How- 
ever, for this behavior to be part of a viable scenario, the dynamics has to be compatible 
with the present-time observational data and it has to admit an early matter-dominated 



phase in order to guarantee structure formation. For the analytic solution (31) and (32) 
these requirements were satisfied. The situation is less clear for the results of the dynamical 
system analysis of the last section. To better understand this point, it is useful to write the 
balances for the components in the form 

p'^ + Pm[l- Tp-^-V"-^ (1 + ry^'] = (58) 
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and 



p', + p,[l + w + Tp'^-V" (1 + r)^+^] = 



(59) 



Of particular interest for an acceptable dynamics is a final attractor (or a stable focus) for 
which both p and r approach their final stationary values from p > Pc and r > Vc respectively. 
This implies the requirements p' < and r' < 0, respectively, during the cosmic evolution. 
The consequences are 

p' <0 =^ r > |w| - 1 = Tc (60) 



and 



r' < 



-fp""-^ < |w|r^-'^(l + r) 



-(s+2) 



(61) 



The condition (61) puts upper limits on the interaction constant 7. For the present epoch 

-{s+2) ^ The right-hand side of the last 



we have tq > |w| — 1 and 'jp^~^ < |u;|rQ~" (1 + tq) 
inequality is of the order of one. The present energy density po is of the order of the 
critical density pcro = 1.88/;.^10~^^gcm"^. Consequently, 7 ^ (lO^^g'^cm'^)™' ^. (Assuming 
tentatively the applicability of the latter inequality also for the analytic solution ( [3l] ) which 
has 171 = 1, this condition reduces to 7 < 1, consistent with the result 7 ^ 0.36 of our data 



analysis.) It follows from (58) and (59) that for any value of 7, considerably smaller than 



this one, the interaction would only provide a very small correction to the dynamics at the 
present epoch. On the other hand, p' < implies 



\w 



i\w\ 



7 



(62) 



This provides us with 



\w\-^'+^U\w\ 



To) 



-(s+2) 



(63) 



and 



(kl-1) 



1— n 



< {iwr'ip^-') . 



(64) 



The inequality (64) relates the interaction strength to deviations from |w| = 1 (cf. Eq. (30) for 
a similar feature of the analytic solutions (31) and (32)). Except for n = 1, any 7P™^^ ^ 1, 



equivalent to a small influence of the interaction on the present cosmological dynamics, 
requires a value w of the EoS parameter close to w = —1. In other words, deviations from 
w = —1 are a measure of the interaction strength. At the same time, a small deviation 
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from w = —1 corresponds to a small, but non-zero limiting value Tc < 1 of the energy- 
density ratio. Since a phantom-type EoS parameter close to -1 is preferred by a number of 
investigations (see, e.g., |13] for a recent analysis), a scenario with finite long-time limits Tc 
and pc does not seems to contradict current observations. 

Let us now consider qualitatively the conditions for the existence of an early matter 
dominated epoch. According to (8), a matter era with p oc requires ^ 1 for a ^ 1. 
For a constant value of the EoS parameter \w\ this is achieved for r ^ 1 at a ^ 1. Since, 



according to (10), the far-future limit of r is smaller then unity (for values of w slightly but 
not substantially smaller than —1), one has r' < as already mentioned, i.e., a decaying 
ratio of the energy densities during the cosmic evolution. Assuming accordingly, that at 



high redshifts r ^ 1 is valid, the balance equation ( 58 ) can be written 



^ ~ _ [1 _ 7r'"+"+'^- V^"'] 

Pm 

Then, under the condition r ^ 1 and with m > 1, the interaction term is negligible for 



[1 _ . (65) 

Pm 



m + n + s < 1 in (65), which corresponds to a matter dominated phase. By inspection one 
realizes that this requirement is met by the first, second and fourth cases of Table [T] and by 
the first, second and third cases of table |Tlj These examples are potential candidates for a 
scenario which correctly reproduces an early matter dominated period but predicts a future 
evolution towards an attractor or a stable focus with a finite, non-vanishing value for the 
ratio of the energy densities. Also the first, second and fourth examples of Table \LU\ have 
m + n + s < 1. For the remaining cases which have m + n + s = 1 there is no dependence on r 
and a matter dominated phase requires 7P^~^ ^ 1 for a ^ 1. Although these considerations 
remain on a heuristic level, its preliminary conclusions do not seem to be inconsistent with 
a scenario that includes standard structure formation. 



VI. DISCUSSION 

Models with an interaction between dark matter and dark energy have received consid- 
erable attention since they provide a framework to address the cosmic coincidence problem. 
We have analyzed here an interacting two-component system of dark matter and dark en- 
ergy with the total energy density p and the ratio of the energy densities of dark matter 
and dark energy r = — as independent variables. We have shown that asymptotically 



19 



finite, positive critical values Vc and pc of these variables require a phantom-type equation 
of state for the dark energy component and an energy transfer from dark energy to dark 
matter, at least in the critical point. These results hold quite generally and do not depend 
on the details of the interaction. We have investigated in some detail interactions of the 
type Q = SH'yp'^'^'^p^plr"' which generalizes a number of models that have been studied so 
far in the literature. For various combinations of m, n and s the resulting dynamics was 
shown to be asymptotically different from that of the ACDM model. As a consequence of 
the interaction, a phantom-type equation of state of the dark energy does not result in a 
big-rip singularity. For the particular interaction Q = 3H'~fp'^/p we found an analytic solu- 
tion that reproduces an early matter-dominated phase and approaches a finite, stationary 
value for the ratio of the energy densities of DM and DE. This solution is consistent with the 
observational data of SNIa from the Union2 data set. It implies a direct relation between the 
interaction strength and deviations from w = —1. The interaction for this case is negligible 
at high redshifts but crucially influences the future dynamics. 

Based on a dynamical system analysis for a broader class of interactions, we identified 
attractors, stable focuses and centers as potential limiting configurations of the cosmological 
dynamics. For attractors and stable focuses (examples are summarized in Tables |l] and |ll| 
respectively,) the asymptotic value gc = — 1 of the deceleration parameter is approached by 



a power of the scale factor. For centers (see the examples in Table III) we find oscillations 
of the deceleration parameter around this critical point, equivalent to a periodic crossing 
of the phantom divide for the EoS of the total cosmic substratum. All these models have 
necessarily m > 1. A qualitative discussion of the conditions under which this class of 
models admits the existence of an early matter dominated phase singles out models with 
■m + n + s < 1. A more quantitative analysis as well as a study of the perturbation dynamics 
of non-linearly interacting models will be the subject of future research. 
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